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Equations of motion, in cylindrical co-ordinates, for the observed rotation of gases within the 
gravitational potential of spiral galaxies have been derived from Carmeli's Cosmological General 
Relativity theory. A TuUy-Fisher type relation results and rotation curves are reproduced without 
the need for non-baryonic halo dark matter. Two acceleration regimes are discovered that are 
separated by a critical acceleration « 4.75 x 10"^" m.s~^ . For accelerations larger than the critical 
value the Newtonian force law applies, but for accelerations less than the critical value the Carmelian 
regime applies. In the Newtonian regime the accelerations fall off as r~^, but in the Carmelian regime 
the accelerations fall off as r~^. This is new physics but is exactly what is suggested by Milgrom's 
phenomenological MOND theory. 
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I. INTRODUCTION 

The rotation curves highlighted by the circular mo- 
tion of stars or more accurately characterised by the 
spectroscopic detection of the motion of neutral hydro- 
gen and other gases in the disk regions of spiral galaxies 
have caused concern for astronomers for many decades. 
Newton's law of gravitation predicts much lower orbital 
speeds than those measured in the disk regions of spiral 
galaxies. 

The most luminous galaxies show slightly declining ro- 
tation curves (orbital speed vs radial position from nu- 
cleus) in the regions outside the star bearing disk, com- 
ing down from a broad maximum in the disk. Inter- 
mediate mass galaxies have mostly nearly flat rotation 
speeds along the disk radius. Lower luminosity galaxies 
usually have monotonically increasing orbital velocities 
across the disk. 

The traditional solution has been to invoke halo 'dark 
matter' 0| that surrounds the galaxy but is transpar- 
ent to all forms of electromagnetic radiation. In fact, 
astronomers have traditionally resorted to 'dark matter' 
whenever known laws of physics were unable to explain 
the observed dynamics. 

In 1983 Milgrom introduced his MOND [H [H [H, 
an empirical approach, which attempts to modify New- 
tonian dynamics in the region of very low acceleration. 
Newton's law describes a force proportional to r~^, where 
r is the radial position from the center of the matter dis- 
tribution, but Milgrom finds that a law fits the data 
very well 

Carmeli 0] approached the problem from a differ- 
ent perspective. He formulated a modification, actually 
an extension of Einstein's general theory, in an expand- 
ing universe taking into account the Hubble expansion as 
a fundamental axiom, which imposes an additional con- 
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straint on the dynamics of particles [3- 

Carmeli believes the usual assumptions in deriving 
Newton's gravitational force law from general relativity 
are insufficient, that gases and stars in the arms of spiral 
galaxies are not immune from Hubble flow. As a conse- 
quence a universal constant oq (in this case, slightly dif- 
ferent to Milgrom's) is introduced as a characteristic ac- 
celeration in the cosmos. Using this theory Carmeli suc- 
cessfully provided a theoretical description of the TuUy- 
Fisher law Q. 

Following Carmeli's lead, Hartnett showed that the 
same line of reasoning leads to plausible galaxy rotation 
curves. The latter used a density model for spiral galax- 
ies, that assumed that most of the mass of the galaxy 
was in the nuclear bulge and that the density in the disk 
region was constant. However by perturbing the den- 
sity one could also get the variation in rotation speeds 
as typically observed. That paper also assumed spheri- 
cal coordinates and a hyperbolic density distribution of 
disk matter, conditions which are not appropriate for 
exponential-density- model galaxies . 

In this paper we take the analysis further, and in a 
more rigorous way, model the gravitational potential and 
the resulting forces determining how test particles move 
in the disks of spiral galaxies using cylindrical coordi- 
nates and an exponential density distribution. Two ac- 
celeration regimes are discovered. In one, normal New- 
tonian gravitation applies. In that regime the effect of 
the Hubble expansion is not observed or is extremely 
weak. It is as if the particles' accelerations are are so 
great that they slip across the expanding space. In the 
other, new physics is needed. There the Carmelian met- 
ric provides it. In this regime the accelerations of parti- 
cles are so weak that their motions are dominated by the 
Hubble expansion and as a result particles move under 
the combined effect of both the Newtonian force and a 
post-Newtonian contribution. 
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II. GRAVITATIONAL POTENTIAL 

In the weak gravitational limit, where Newtonian grav- 
itation applies, it is sufficient to assume the Carmeli 
metric with non-zero elements 500 = f + (744 = 

1 + gkk = —1, {k = I, 2, 3) in the lowest approxi- 

mations in both I/c and I/t. Here a new constant, called 
the Hubble-Carmeli constant, is introduced r « I/Hq. 
See @ for details. The potential functions cj) and are 
determined by Einstein's field equations and from their 
respective Poisson equations, 

V20=:47rGpe//, 
47rGpe// 



(la) 
(lb) 



where p is the mass density and ao a universal charac- 
teristic acceleration gq — c/t. As usual c is the speed of 
light in vacuo. 

In Carmelian theory p^f / = p — Pc is used instead of 
matter density p. The parameter p^ is the critical density 
of the universe given by pc = S/SttGt^ « 10~^^ g.cm~^. 
However in a galaxy because p ^ Pc, Pc can be neglected 
in this paper. 

A comparison of and ip in Hla|l and ljlb|l leads to 
■0 = (t>/cLo within an arbitrary additive constant. Since 
both potentials are defined with respect to the same co- 
ordinate system, in reality, we only need deal with one 
potential function, the gravitational potential, (j). 

In cylindrical coordinates (r, d-^_z) the potential </) that 
satisfies ()la|l can be found from [lj]i 



r) = ~27rG / Jo{kr)dk 

Jq Jo 



p{r')Joikr')r'dr' (2) 



where Jo{kr) is the zeroth order Bessel function function 
and k is the z coordinate scale factor {k = 1/b). It is 
also assumed that the density function can be modeled 
as a delta function of the vertical coordinate z. There- 
fore the density p{r, z) — p{r)p{z) — p(r)S{z) with no 9 
dependence. To correctly model the effect of the spiral 
arms a 9 dependence may be needed, but for our model 
it is assumed independent. The requirement on the z de- 
pendence is satisfied with density functions of the form 



p(z) = -isec/i(^)2 or e"l"l/^ 



t <=" 

Here, provided the scale length b is much smaller than the 
limit of the actual matter distribution in the z direction 
then the integral over all z yields a contribution to the 
mass of unity. This is the thin disk approximation which 
seems to be fairly applicable over both disk and galactic 
bulge. 

The integral over dk in ((SJ is the surface density which 
may be calculated once the form of the density p is 
known. Following from observation we choose an expo- 
nential function of the form 

M 
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FIG. 1: TuUy-Fisher law plotted on logarithmic axes. Curve 
1 (solid line) represents the fourth order dependence of the 
rotational speeds of tracer gases in galaxies determined from 
the Carmelian equation 1121 . The masses are expressed in 
solar mass units of Mq — 2 x lO^^ftp. Curve 2 (broken line) 
represents the straight line = 2 x 10^ + 0.064M (Mass) 



for the radial dependence, where a is a radial scale length 
and M is the mass of the galaxy. 



III. EQUATIONS OF MOTION 

The Hubble law describes the expansion of the cosmos 
and the matter embedded in it. Therefore the line ele- 
ment for any two points in this 'new' space-time-velocity 
is ds^ — gmc^dt^ + gkk{dx^Y + gur'^dv'^ = 0. Here 
k — 1,2,3. The relative separation in 3 spatial coor- 
dinates — (x^)^ -I- (a;^)^ -I- (x^)^ and the relative ve- 
locity between points connected by ds is v. The Hubble- 
Carmeli constant, r, is a constant for all observers at the 
same epoch, therefore may be regarded as a constant on 
the scale of any measurements. 

The equations of motion (B.62a and B.63a from |^) to 
lowest approximation in 1/c are reproduced here, 



d'^x^ 
~dP' 



1 dcj) 



(5) 



This is the usual looking geodesic equation derived from 
general relativity but now in 5 dimensions. And the sec- 
ond is a new 'phase space' equation derived from the 
Carmeli theory, 



d X 
dv"^ 



1 av 

'2dx^' 



(6) 
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FIG. 2: (a) Above: The rotational speeds of tracer gases in 
NGC 3198 (SBc barbed spiral) (circles - curve 1). Theoretical 
curve fits from the Carmelian equation II15II (broken curve 2) 
and from the Newtonian equation (|SJ (curve 3) (b) Below: 
The critical acceleration |ao (curve 1). The rotational ac- 
celerations determined from the Carmelian (curve 2) and the 
Newtonian (curve 3) equations with their respective values of 
a and M 
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FIG. 3: (a) Above: The rotational speeds of tracer gases in 
NGC 0598 (Sc spiral) (circles - curve 1). Theoretical curve 
fits from the Carmelian equation 1)15^ (curve 2) and from the 
Newtonian equation Q (curve 3). (b) Below: The critical 
acceleration |ao (curve 1). The corresponding rotational ac- 
celerations determined from the Carmelian (curve 2) and the 
Newtonian (curve 3) equations 



Newtonian 



It follows from (|SJ|, and (|2J), and the usual form of 
the circular motion equation 



V 

r 



dr ' 



that 



2 GMv^ r. 

where G denotes the gravitational constant and 



(7) 



(8) 



where / and K are standard zeroth and first order Bessel 
functions. 

Equation ^ is the usual Newtonian result for the 
speed of circular motion in a cylindrical gravitational po- 
tential. This equation has been plotted in curve 3 of figs 
|2Ia) - llUf a) as a function of radial position from the cen- 
ter of a galaxy in kiloparsecs (kpc) where kpc w 3 x 10^® 
m. Best fits were determined with M and a as free pa- 
rameters. Through out this paper M is expressed in solar 
mass units M© « 2 x 10'^" kg. 



B. Carmelian 

Using ip — (p/ttQ in © results in a new equation 

dr 

V = ao 



(10) 



where we have integrated and solved for t; as a function 
of r. Using the potential cj), determined from Q and 
in H10|) . results in 



2 r3/2 

V = -ao r—— -, 

3 VGM 



(11) 



which describes the expansion of space within a galaxy. 

It must be realised that the only direction in the cylin- 
drical coordinates of a galaxy that is free to expand in 
the Hubble flow is the azimuthal Therefore Hll|) de- 
scribes the velocity component in that direction. Carmeli 
applied this line of reasoning 0, |^ . 

To establish the combined result of the two equations 
of motion © and ((TT|l . the simultaneous speed of test 
particles must be determined by the elimination of r be- 
tween the two equations. The physical meaning can be 
understood in terms of particles that simultaneously sat- 
isfy both © and © . 
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FIG. 4; (a) Above: The rotational speeds of tracer gases in 
NGC 2903 (Sc spiral) (circles - curve 1). Theoretical curve 
fits from the Carmelian equation l|15^ (curve 2) and from the 
Newtonian equation Q (curve 3) (b) Below: The critical ac- 
celeration |ao (curve 1). The corresponding rotational ac- 
celerations determined from the Carmelian (curve 2) and the 
Newtonian (curve 3) equations 




Radial position [kpc] 

FIG. 5: (a) Above: The rotational speeds of tracer gases in 
NGC 7331 (Sbc spiral) (circles - curve 1). Theoretical curve 
fits from the Carmelian equation 1)15^ (curve 2) and from the 
Newtonian equation Q (curve 3) (b) Below: The critical ac- 
celeration |ao (curve 1). The corresponding rotational ac- 
celerations determined from the Carmelian (curve 2) and the 
Newtonian (curve 3) equations 



The usual Newtonian expression ijSJ) describes motion 
under the central potential but assumes that spatial co- 
ordinates are fixed. Whereas the new equation © de- 
scribes the expansion of space itself within a galaxy. 
Therefore we must find the combined (simultaneous) ef- 
fect of © and Hll|) . The result is a post-Newtonian equa- 
tion, 



„2/3 



(GM) 



5/3 



n, 



(12) 



derived from (|HJ) where the following substitution 

has been made from l|ll(l . The resulting equation, here- 
after referred to as Carmelian, cannot be solved analyti- 
cally. However using the Mathematica software package 
(|12|l can be solved numerically. 

The result is plotted in curve 1 of fig. ^ where it has 
been assumed that a = 1 kpc and it is compared with 
the straight line v"^ = 2 x 10^ + 0.064M (curve 2). For 
large M the small offset can be neglected. This result 
indicates that the fourth order dependence on rotational 



speed (v) is directly proportional to mass (Af ) for large 
masses. 

Assuming that the masses, of the galaxies studied, are 
directly proportional to their luminosity, this dependence 
then becomes the Tully-Fisher relationship. This extends 
the work of Carmeli and derives the underlying the- 
oretical framework upon which the Tully-Fisher law is 
founded. 



C. Rotation curves 

In 0,|EI1|> using spherical co-ordinates, it was found 
that in the limit of large r and where all the matter was 
interior to the position of a test particle, such a particle 
is also subject to an additional circular motion described 
by (|ll|l . Apparently this is the result of the expansion of 
space itself within the galaxy but in an azimuthal direc- 
tion to the usual center of coordinates of the galaxy. In 
this paper also the same result (|11|) was obtained but in 
this case derived using cylindrical coordinates. 

Carmeli 3j determined a Tully-Fisher type relation 
using the Newtonian circular velocity equation expressed 
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FIG. 6: (a) Above: The rotational speeds in NGC 2841 (Sb 
spiral) (circles - curve 1). Theoretical curve fits from the 
Carmelian equation H15|l (curve 2) and from the Newtonian 
equation ||HJ (curve 3) (b) Below: The critical acceleration 
|ao (curve 1). The corresponding rotational accelerations 
determined from the Carmelian (curve 2) and the Newtonian 
(curve 3) equations 



FIG. 7: (a) Above: The rotational speeds in IC 0342 (Sc 
spiral) ((circles - curve 1). Theoretical curve fits from the 
Carmelian equation l|15|l (curve 2) and from the Newtonian 
equation ||HJ (curve 3) (b) Below: The critical acceleration 
|ao (curve 1). The corresponding rotational accelerations 
determined from the Carmelian (curve 2) and the Newtonian 
(curve 3) equations 



in spherical coordinates, 



GM 



where it is assumed that test particles orbit at radius 
r outside of a fixed mass M. Then by eliminating r 
between (|13|) and Hll|l we get the result. This is achieved 
by taking the 3/2 power of (|13() and multiplying it with 
(fTT|l yielding 



= GM-ao. 



which is the radial position (or r) dependent part of H15|l . 
Hence (jl5|l then recovers the form of the TuUy-Fisher 
(-23~) relation ffT^ . 

By taking the 4th root of (|15|l we get an expression 
for the circular velocity of test particles as a function of 
their radial position r. That result has been plotted in 
curve 2 of figs [21 a) - llUf a) for each galaxy with a and M 
determined as fit parameters. The resulting curves have 
the characteristic flat shape for large radial position r. 
At small values of r the rotation speeds determined from 
the Newtonian equation (curve 3) dominate as seen in 
(14) figsHa) -Cnia). 



So by applying the same approach with (jH)) (raising it 
to the 3/2 power) and multiplying it with we can de- 
rive an equation describing the rotation curves in galax- 
ies. The result is 



= GM-ao 



r \9/2 
2a/ 



n3/2 



(15) 



remembering 11 is a function of r/2a. It is is easily con- 
firmed that as r — > cxD, 



(2a) 



9/2 



n3/2 



IV. ACCELERATIONS 

The acceleration |ao in (|ll|l can be considered to be 
a critical acceleration. Therefore when we compare the 
accelerations derived from the Newtonian equation © 
and the Carmelian equation (|15l) with this critical accel- 
eration we notice two regimes develop. See figs EJb) - 
CUb). 

For example, fig ^b) is very instructive. There the 
straight line (curve 1) is the critical acceleration |oo, 
curve 2 represents the acceleration derived from the 
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Carmelian equation H15() and curve 3 represents the accel- 
eration derived from the Newtonian equation © . For the 
values, determined from the fits, of the mass (M) and the 
radial scale length (a), which determine how the matter 
density varies as a function of radial position r, the curves 
2 and 3 cross each other very close to the critical acceler- 
ation. The significance is that for accelerations less than 
the critical acceleration the Carmelian force applies and 
for accelerations greater than the critical acceleration the 
Newtonian force applies. Note also that the Newtonian 
curve 3 has a dependence and the Carmelian curve 2 
has a dependence above 10 kpc. When curves of the 
form were fitted to the functions used in figEIb) be- 
tween 15 and 20 kpc, the coefHcients x were determined 
to be x = 2.003 and x = 1.025 respectively. 

From ||SJ) the gravitational acceleration [v^ /r) can be 
calculated in the limit of r — > c», outside most of the 
matter of the galaxy. As expected for the Newtonian 
model it tends to GM/r^ . And similarly from H15|l the 
gravitational acceleration (u^/r) can also be calculated 
in the limit of large r, for the Carmelian model. In 
this case it is evident from (|14|l that it must tend to 

^JcM^oo/r. In this regime the accelerations are very 
weak. This is very significant as alternative theories of 
gravity have been suggested (for example, [Tol UlL ll^l 
where the force of gravity falls away as for small ac- 
celerations. However for small r, that is, where r — > 0, 
close to the origin of the central gravitational potential, 
the effect of the Carmelian force law becomes extremely 
small and is many orders of magnitude smaller than that 
for the Newtonian force law. 
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FIG. 8: (a) Above: The rotational speeds of tracer gases in 
NGC 1097 (SBb barbed spiral) (circles - curve 1). Theoreti- 
cal curve fits from the Carmelian equation 1151 (curve 2) and 
from the Newtonian equation (|HJ (curve 3) (b) Below: The 
critical acceleration |ao (curve 1). The corresponding rota- 
tional accelerations determined from the Carmelian (curve 2) 
and the Newtonian (curve 3) equations 



A. Extragalactic spirals 



V. SAMPLE OF GALAXY ROTATION CURVES 

A sample of 9 galaxy fits are shown in figs [5] -El The 
top (a) figures show the rotation curve fits and the bot- 
tom (b) figures show the resulting acceleration regimes. 
In each figure, the measured rotational speeds of tracer 
gases in the chosen spiral galaxy is shown func- 
tion of radial position (curve 1). Measured data are taken 
from 13]. Theoretical curves from the Carmelian equa- 
tion (|15|) (curve 2) and from the Newtonian equation © 
(curve 3) are fitted over the range of r which best fit the 
data by allowing a and AI to be free parameters. The 
accelerations in the bottom (b) figures are derived from 
the Carmelian (curve 2) and the Newtonian (curve 3) 
equations respectively with values of a and M derived 
from the fits in the (a) figures. These are compared with 
(curve 1) the critical acceleration |ao « 4.75 x 10~^° 
m.s~^ determined elsewhere from r w 4.21 x lO^^s, which 
is the reciprocal of the Hubble parameter at zero distance 
/i« 73.270 km s-^Mpc-^. 



In this section, I discuss individual galaxy curve fits, 
starting with fig El showing the barbed spiral NGC 3198. 
In each case, because of the possibility of different accel- 
eration regimes, both Carmelian and Newtonian curve fit 
were attempted. In fig Ofa) the Carmelian fit is shown 
by the broken curve 2 to be the only good fit. The scale 
radius a = 1.85 kpc and M = 0.984 x W^Mq deter- 
mined from the fit. Curve 3 shows the best Newtonian 
fit with a = 2.99 kpc and M = 4.2 x IO^^Mq but it 
doesn't fit well at high values of r. The Newtonian fit 
results in a mass at least 4 times greater than that from 
the Carmelian fit. The scale radius determined from lu- 
minous matter for this galaxy is a = 2.5 kpc which is 
closer to the Carmelian curve determination. 

In the figl3b) curve 2 shows the acceleration using the 
Carmelian determined values of the scale radius a and 
mass M. Curve 3 shows the acceleration for Newtonian 
fit determined values. Clearly curve 2 is dominant and 
is always less than the critical acceleration |ao. In this 
model, when the accelerations are less than the critical 
value the Carmelain force applies. 

Next the data for NGC 0598, a Sc spiral galaxy, shown 
in figOIa), fits both a Carmelian and a Newtonian curve. 
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FIG. 9; (a) Above: The rotational speeds of tracer gases in 
NGC 2590 (Sb spiral) (circles - curve 1). Theoretical curve 
fits from the Carmelian equation l|15^ (curve 2) and from the 
Newtonian equation Q (curve 3) (b) Below: The critical ac- 
celeration |ao (curve 1). The corresponding rotational ac- 
celerations determined from the Carmelian (curve 2) and the 
Newtonian (curve 3) equations 



From the Carmelian model the scale radius a = 1.85 kpc 
and M = 0.13 X 10^°Mq determined from the fit. Clearly 
the Carmelian curve 2 is the better fit over the Newtonian 
curve 3 for the following reasons. Firstly the data (cir- 
cles - curve 1) continue to rise or at least are not falling 
at the extremity of the available measured range. The 
Newtonian curve indicates it should fall. Secondly from 
figEIb) the accelerations are much less than the critical 
value |ao and hence in this regime the Carmelian force 
law applies. 

From the Newtonian model a scale radius of a = 
2.22 kpc and M = 1.42 x 10^°Mq were determined but 
the fit doesn't conform to the model. Nevertheless the 
Newtonian fit results in a mass 10 times greater than the 
fit for the Carmelian model. 

Figs 0] and |S1 show both Newtonian and Carmelian 
models fit the rotation curve data for the galaxies NGC 
2903 and NGC 7331. In the high acceleration regime 
a Newtonian fit is applicable and when the acceleration 
drops below |ao the Carmelian applies. From the Table 
the best fit determined values of a and M for each galaxy 
are listed. They are compared with published values of 
a and masses determined from different methods. 

In these and all subsequent top (a) figures, the fits 
for curves 2 and 3, respectively, apply only for acceler- 
ations less than and greater than the critical accelera- 



tion. In the rotation curve fits, in the top (a) figures, 
the circular velocities are not added but each apply over 
their respective regimes. This means the masses deter- 
mined from the Newtonian fits must be less than those 
from the Carmelian fits because the Newtonian deter- 
mined mass must lay within the radius R where the two 
curves cross. The Newtonian fits occur in the stronger 
acceleration regimes (> |ao) close to the galactic center 
as indicated by all of the bottom (b) figures. 

The Table lists the radial position (i?) where the New- 
tonian and Carmelian regimes meet, the mass (M^) out 
to a radius r = R and the total mass M shown in Col. 
8 (Carmeli M from fit), which is the Carmelian model 
determined mass. As expected the values of Mr < M 
determined from the Newtonian regime. 

Figs El and d respectively, show NGC 2841 and IC 
0342, which have been modeled with two central com- 
ponents with different scale radii. Both seem to have a 
dense mass concentration toward their centers. The scale 
radii for these inner most concentrations are a = 0.09 kpc 
and a = 0.05 kpc for NGC 2841 and IC 0342 respectively. 
The combined rotation curve for IC 0342 is a much better 
fit to the data using the two-acceleration-regime model 
than that of the former, NGC 2841. Nevertheless the 
theory works well in these type of galaxies also. 

Figs El and El respectively, show NGC 1097 and NGC 
2590 modeled with only one central dense Newtonian 
component and they fit the model very well. Deviations 
in most cases I believe can be attributed to the fact that 
no azimuthal dependence has been added to the model, 
nor does it allow for anything but constant scale radii 
over regions of the galaxy where the fits apply. That is 
clearly unphysical but seems to be a reasonable approx- 
imation. 

The Table also shows, in Col. 9, the mass (Afio) de- 
termined at r = 10 kpc using the Newtonian formula 



(16) 



assuming that most of the mass is internal to r = 10 kpc 
and the measured rotation speeds. This calculation is 
compared with the mass (Mio), in Col. 10, derived from 
the Carmelian equation with r = 10 kpc. The ratios of 
these masses are shown in Col. 11. It indicates that 
using normal Newtonian/Keplerian dynamics seriously 
over estimates galaxy masses by between 2 and 7 times. 
These values are consistent with the needed mass levels of 
halo dark matter to achieve the correct rotation curves. 



B. The Galaxy 

Considering the (Milky Way) Galaxy the same analy- 
sis has been applied to the data from 0. See the rota- 
tion curve in fig llOf a). Other observers don't record the 
central peak indicative of a large central mass concen- 
tration. However the compact radio source Sagittarius 
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Table: Important galaxy data where valid curve fits are found 









Scale radius 


Newton 


Carmeli 


Newton 


Carmeli 


Newton 


Carmeli 




Radius 




Figure 


Galaxy 


Tvpe 


a 


a 


a 


M 


M 


Mio 


Mio 


Ratio 


R 


Mr 




Name 




published 


from fit 


from fit 


from fit 


from fit 












2 


NGC 3198 


SBc 


2.5 




1.85 




0.984 


6.554 


0.956 


6.85 






3 


NGC 0598 


Sc 


2.7 




1.01 




0.13 












4 


NGC 2903 


Sc 


1.9 


0.31 


0.98 


0.54 


2.12 


8.810 


2.120 


4.16 


1.39 


0.54 


5 


NGC 7331 


Sbc 


4.7 


1.15 


1.2 


4.0 


4.45 


11.496 


4.440 


2.59 


3.38 


3.16 


6 


NGC 2841 


Sb 


2.3 


2.26 


2.30 


8.12 


9.00 


17.721 


8.377 


2.11 


2.11 


5.09 


7 


IC 0342 


Sc 




0.74 


1.05 


1.28 


1.81 


8.200 


1.809 


4.53 


2.08 


1.00 


8 


NGC 1097 


SBb 




0.62 


2.12 


4.80 


9.74 


22.997 


9.242 


2.49 


3.94 


4.74 


9 


NGC 2590 


Sb 


2.1 


0.5 


2.73 


2.90 


6.96 


17.213 


6.127 


2.81 


3.73 


2.89 


10 


MW Galaxy 


Sb 




0.12 


1.09 


0.45 


2.31 


9.302 


2.308 


4.03 


1.22 


0.45 



Cols (4) -(6): a in kpc. Cols (7)-(10), (13): M in [10^° M©]. Cols (9) & (10): Mass (Mio) calculated at r = 10 kpc. Col.(ll): 
Mass from Col. (9) divided by Col. (10). Cols (12): R in kpc. Col. (13): Mass {Mr) calculated ai r = R. 
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range 3 to 10 kpc is an excellent fit. The acceleration 
regimes, consistent with the curve fits (curves 2 and 3) 
in fig llOf a). are shown in fig llOf b). Except for the in- 
tervening region the fits agree well with the theory. The 
discrepancy could be due to an unmodeled higher concen- 
tration of mass in the central bulge region. Deviations in 
the spiral arms from a smooth exponential density decay 
are consistent with the oscillations about curve 2 in fig 
da). 

The distance of our solar system from the Galactic 
center has been estimated at 9.95 kpc |^ and more re- 
cently as 7.94 ± 0.42 kpc 0- Likewise the enclosed mass 
may be calculated using both IjlGI) and (|15|l . When one 
uses the normal Newtonian/Keplerian calculation 1)16(1 
and V = 200.78 A:to.s~^ the speed of the solar system 
orbiting the Galactic center, it results in an estimate 
of the enclosed mass of Mio — 9.3 x lO^^M© as com- 
pared with Mio = 2.3 X lO^^M© from the Carmelian 
calculation H15|l . which is 4 times smaller. For a distance 
7.94 kpc the Keplerian calculation of the enclosed mass 
is Mio = 7.6 X 1O^°M0, which is 3.3 times greater than 
the Carmelian calculation for the enclosed mass at that 
distance. 



FIG. 10: (a) Above: The rotational speeds of tracer gases 
in the Galaxy (Milky Way Sb spiral). Measured data from 
fl^ shown in curve 1. Theoretical curve from the Carmelian 
equation (|15|l (curve 2) and from the Newtonian equation 
^ (curve 3) (b) Below: The critical acceleration |ao (curve 
1). The corresponding rotational accelerations determined 
from the Carmelian (curve 2) and the Newtonian (curve 3) 
equations respectively 



A* at the Galactic center is widely believed to be asso- 
ciated with the supermassive black hole with a mass of 
about 3.59 ± 0.59 x IO^Mq. Then there is the vast 
concentration of matter in the Galactic bulge. From the 
Newtonian curve fit we'd expect that within 1.2 kpc of 
the center there is a mass of about M = 4.5 x lO^M©. 
The Carmelian curve fit (curve 2 in fig llOf a)) over the 



VI. CONCLUSION 

Carmeli's Cosmological General Relativity theory pro- 
vides a solution to the rotation curve anomaly in the 
outer regions of spiral galaxies. Equations of motion have 
been derived from Carmeli's metric assuming a gravi- 
tational potential in cylindrical coordinates. A TuUy- 
Fisher type relation results and the rotation curves for 
spiral galaxies are accurately reproduced without the 
need for non-baryonic halo dark matter. 

Two acceleration regimes are discovered that are sep- 
arated by a critical acceleration |ao. For accelerations 
larger than the critical value the Newtonian force law 
applies, but for accelerations less than the critical value 
the Carmelian regime applies. In the Newtonian regime 
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the accelerations fall off as as expected, but in 
the Carmelian regime the accelerations fall off as r~^. 
This is new physics but is exactly what is suggested by 
Milgrom's phenomenological MOND theory. However in 
this case a theoretical basis is found, whereas until now, 
no theory has been found for Milgrom's MOND. This 
theory also provides an understanding of the connection 
between the two regimes. 
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